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■ ABSTRACT 

O '• 
O ; 

I BIMA observations of the Orion nebula discovered a giant flare from a young 

^! star previously undetected at millimeter wavelengths. The star briefly became 

^ ! the brightest compact object in the nebula at 86 GHz. Its flux density increased 

1^ ! by more than a factor of 5 on a timescale of hours, to a peak of 160 mJy. This 

^ ! is one of the most luminous stellar radio flares ever observed. Remarkably, the 

^ i Chandra X-ray observatory was in the midst of a deep integration of the Orion 

^ i nebula at the time of the BIMA discovery; the source's X-ray flux increased by 

i a factor of 10 approximately 2 days before the radio detection. Follow-up radio 

^ i observations with the VLA and BIMA showed that the source decayed on a 

O ! timescale of days, then flared again several times over the next 70 days, although 

2 I never as brightly as during the discovery. Circular polarization was detected 

^ I at 15, 22, and 43 GHz, indicating that the emission mechanism was cyclotron. 

VLBA observations 9 days after the initial flare yield a brightness temperature 

O I Tfe > 5 X 10^ K at 15 GHz. Infrared spectroscopy indicates the source is a K5V 

I star with faint Br 7 emission, suggesting that it is a weak-line T Tauri object. 

I Zeeman splitting measurements in the infrared spectrum find B ~ 2.6 ± 1.0 kG. 

I The flare is an extreme example of magnetic activity associated with a young 

^ I stellar object. These data suggest that short observations obtained with ALMA 

I will uncover hundreds of flaring young stellar objects in the Orion region. 

Subject headings: open clusters and associations: individual (Orion Nebula Clus- 
ter) — stars: flare — stars: formation — stars: magnetic fields 
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1. Introduction 

The transient radio sky is mostly unexplored, especially at millimeter wavelengths. 
The few systematic radio searches for astronomical transients have been done at longer 
wavelengths, cover a small portion of the sky, or have low sensitivity (e.g., Langston et 
al. 2000, Hyman et al. 2002, Carilli et al. 2003), principally due to the lengthy observing 
necessary to conduct transient surveys with even the most sensitive of existing instruments. 
Many of the known radio transients, in fact, were first discovered at other wavelengths (e.g.. 
Frail et al. 2003). Therefore, the population of transient millimeter wavelength sources is 
currently only poorly constrained. With the advent of new, more sensitive facilities (such as 
ALMA) , we expect that discoveries of radio transients will be a common occurrence. 

Radio stars are among the best studied transient sources at centimeter wavelengths 
(Giidel 2002). A wide range of stellar types have been shown to emit in the radio: brown 
dwarfs, T Tauri stars, dwarf stars, and evolved giants. The radio properties of young stellar 
objects (YSOs) have been studied in a variety of environments, particularly in the Orion 
nebula (Garay, Moran & Reid 1987; Churchwell et al. 1987; Felli et al. 1993a, 1993b). These 
surveys show the existence of sub-arcsecond radio sources, many of which are variable and 
coincident with stars. Detection of radio circular polarization in YSOs has strengthened 
the case that flares from these stars originate in massive plasma ejections related to coronal 
magnetic field activity, qualitatively similar to those seen in the Sun (Feigelson, Carkman 
& Wilking 1998, White, Pallavicini & Kundu 1992). Although deep X-ray and infrared 
surveys of the Orion nebiila have identified many variable sources (Feigelson et al. 2002; 
Hillenbrand & Carpenter 2000; Carpenter, Hillenbrand & Skrutskie 2001), the origins of 
magnetic activity in YSOs are still poorly understood (Feigelson et al. 2003). 

We report here the serendipitous discovery of a millimeter-wave transient source in 
the Orion nebula cluster (Figure 1). The Berkeley-IUinois-Maryland Association Array 
(BIMA) detected the source at 86 GHz on 2003 January 20, with a flux density that varied 
substantially on a timescale of hours (Bower, Plambeck & Bolatto 2003). For a time the 
object became the brightest source in the cluster, with a flux density greater than those of 
the Becklin-Neugebauer Object (BN) and the massive young star IRc2. This detection was 
possible only because BIMA was in its highest resolution configuration, which provides a 0'/5 
synthesized beam and filters out the much stronger extended emission from dust, molecular 
lines, and the Orion HE region. 

The source position, measured to an accuracy of ±0.1" from the BIMA image, was 
found to be coincident with a variable centimeter-wave radio source (GMR-A; Garay et al. 
1987, Felli et al. 1993a,b), a bright near infrared point source (number 573; Hillenbrand & 
Carpenter 2000), and a variable X-ray source (number 297; Feigelson et al. 2002). These 
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previous observations suggested that GMR-A is a YSO deeply obscured by the molecular 
cloud. New observations presented here support and expand this conclusion. 

Follow-up observations with BIMA and the Very Large Array (VLA) showed that the 
source decayed within days of its outburst. GMR-A flared again on several occasions over the 
70 days following the first fiare, although never to the discovery intensity. The Nobeyama Ra- 
dio Observatory (NRO) detection of GMR-A at 100 and 150 GHz a few days after discovery 
(Nakanishi et al. 2003) found an inverted spectrum, indicating that the synchrotron self- 
absorption frequency of the plasma was greater than 100 GHz. Circularly polarized emission 
was detected with the VLA at several epochs, indicating that the radio emission was caused 
by cyclotron radiation from mildly rclativistic electrons in a strong magnetic field. Very 
Large Baseline Array (VLB A) observations within 10 days of the discovery showed that a 
significant fraction of the fiux density of GMR-A was unresolved on the milliarcsecond scale. 

By an extraordinary coincidence, the Chandra X-ray Observatory was in the midst of 
a very deep integration on the Orion nebula at the time of the BIMA discovery (Getman 
et al. 2003). The X-ray fiux from GMR-A increased by a factor of ~ 10 approximately 2 
days prior to the BIMA observation. To our knowledge, this is just the second simultaneous 
radio/mm and X-ray observation of a YSO (Feigelson et al. 1994) and the only simultaneous 
radio/mm and X-ray observation of a flaring YSO. 

Infrared photometry at the Keck 1 10-m and the Cerro Tololo Inter american Observatory 
(CTIO) 4-m telescopes within days of discovery showed no variabihty with respect to the 
historical fiux. Infrared spectroscopy with the Keck II 10-m telescope allowed a robust 
identification of the infrared source as a star of spectral type K5V, with a fast rotation, a 
measurable Zeeman effect indicative of magnetic activity, and a weak Brackett 7 emission 
line. All evidence points to a strong magnetic outburst from a highly obscured young stellar 
object, probably a weak-line T Tauri star. 

2. Observations and Data Analysis 

This paper presents observations obtained at BIMA, the VLA, the VLBA, Keck, and the 
CTIO. Plots of the light curve including dates for infrared photometry, infrared spectroscopy 
and X-ray observations are shown in Figures 2, 3 and 4. In addition to the observations 
described below, we have included the NRO observations reported by Nakanishi et al. (2003) 
and the Chandra X.-rety observations reported by Getman et al. (2003). Positions determined 
from various methods are listed in Table 1. Flux densities from BIMA, VLA, and near 
infrared observations are fisted in Tables 2, 3, and 4. We discuss each data set in detail 
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below. 



2.1. BIMA Observations 

Continuum observations of the Orion nebula at 86 GHz were made with the BIMA array 
on 2002 December 23 and 2003 January 20 as part of a long term project to track the proper 
motion of the Beckhn Neugebauer Object relative to IRc2 (Plambeck et al. 1995). BIMA 
was in its A-configuration (Welch et al. 1996), which provides telescope spacings of up to 
1.9 km; the synthesized beam on Orion was 0'.'9 x O'.'S. The intense v=l J— 2 — > 1 SiO maser 
associated with IRc2 was used as a phase reference for self-calibration. 

Figure 1 shows the discovery image, from the 2003 January 20 data set. The map covers 
a 2^5 X 2(5 region and has not been corrected for attenuation by the 2(2 BIMA primary beam. 
The flaring object is conspicuous in the upper right quadrant. Note that emission from dust, 
molecular lines, and the Orion HII region is almost completely resolved out in this image. 
There is no sign of the flare source in the 2002 December 23 image, to a 3a upper limit of 
6 mJy (after correction for primary beam attenuation). BIMA observed this fleld also on 
2003 January 11 and 14, during very poor weather. Maps generated from these data show 
BN and IRc2 but no sign of the flare source, to levels of 12 and 15 mJy, respectively. 

Following the preset schedule, BIMA was reconflgured from the A-array into the B- 
array on 2003 January 21, so all follow-up observations of the flare were obtained with a 
5'/5 X 2'/5 synthesized beam. To exclude extended emission, only basehnes longer than 20 kX 
were used to make these maps. Fortunately, the flare source is far enough from the Orion- 
KL nebula that confusion from dust or spectral line emission is unimportant at this angular 
resolution. Brief follow-up observations were conducted on 2003 January 23, 24, and 30 and 
2003 February 4 and 6. GMR-A was detected only on 2003 February 6, with a flux density 
of 11 ± 3 mJy. 

2.2. VLA Observations 

VLA observations were carried out on 19 occasions beginning 2003 January 22 (two 
days after the discovery observations) and ending 2003 March 29. Before 17 February 2003 
intervals between observations ranged from one day to 5 days. Observations had durations 
ranging from 1 to 3.5 hours. They were conducted in the standard continuum mode with 
two polarizations and two IF bands of 50 MHz. Observing frequencies of 8.4, 15, 22 and 43 
GHz were used, although not all bands were observed in each epoch. 
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The data were analyzed in AIPS. The absolute flux density scale was set by observations 
of J0713+438. Following the VLA calibrator manual, we assume for J0713+438 5*8.4 = 
1.14 Jy, 5*15 = 0.73 Jy, 5*22 = 0.55 Jy, and 5*43 = 0.29 Jy. These fluxes agree within 
10% of measurements from 27 December 2002 and 8 February 2003 at 8.4, 22 and 43 GHz 
(http://www.aoc.nrao.edu/~smyers/calibration/). Antenna pointing was stabilized using 
X-band reference pointing on J0607-085. Phase, amplitude, and leakage term calibration 
were performed with observations of J0607-085. Since these observations are brief, the 
leakage term calibration cannot be regarded as providing accuracy in the linear polarization 
better than 1%. Instrumental errors in the fractional circular polarization are < 1% (Bower 
et al. 2002). 

The VLA was in the DnC and D configuration during these observations, providing a 
typical resolution of 3" x 1" at 22 GHz. Because of the substantial extended flux associated 
with the Orion nebula, only visibilities on baselines longer than 20 kX were used. The typical 
image rms noise in a 1-hour long observation was 10 mJy, 3 mJy, 1 mJy and 1 mJy at 8.4, 15, 
22 and 43 GHz, respectively. Flux densities for GMR-A were measured by fitting for a point 
source at the phase center. GMR-A was never detected at 8.4 GHz because confusion from 
extended emission in the 10" synthesized beam dominated the source. We also measured 
the flux density and the position of BN at 22 GHz at each epoch. The stability of these 
quantities gives us confidence in the accuracy of our measurements for GMR-A. We find a 
mean flux density for the BN object at 22 GHz of 12.9 ± 0.1 mJy. 



2.3. VLBA Observations 

Observations with the VLBA were carried out on 2003 January 24 and 29 for 2 and 6 
hours, respectively. Observing frequencies were 15 and 22 GHz in both experiments. The 
observations were performed in a phase-referenced mode that allows detection of a weak 
source. The reference source used was J0541-0541, which is 1.6 degrees away from GMR-A. 
We assumed a position a =05:41:38.0833740 and S =-05:41:49.428471 (J2000) for J0541- 
0541. A cycle time (that includes a cahbrator observation, slew to the source, a source 
observation and slew back to the calibrator) of 40 seconds was used in the first epoch. Cycle 
times of 60 and 120 seconds were used at 15 and 22 GHz in the second epoch, respectively. 
Fringe solutions were determined for most antennas in both experiments using the calibrator 
J0541-0541. One iteration of self-calibration amplitude solutions were also computed for 
J0541-0541, which substantially improved the image quality. These solutions were applied 
to GMR-A. 

A large field was imaged and the source was clearly identified in the second epoch at 
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15 and 22 GHz ~ 40 mas from the field center (Figure 5). Exact coordinates are given in 
Table 1. The reported errors in the position are ~ 20 /xas in right ascension and ~ 40 /xas 
in declination. These are statistical errors. Following Reid et al. (1999), we estimate that 
systematic errors due to miscalibration of the opacity are ~ 100 /xas. 

During the second VLBA epoch, the integrated flux densities measured for GMR-A were 
11.2 ± 0.8 mJy and 16.7 ± 1.3 mJy at 15 and 22 GHz, respectively. These are substantially 
less than the VLA fluxes measured 1 day later, of 32.4 ±1.7 mJy and 30.1 ± 0.8 mJy at 
15 and 22 GHz. The source appears point-like in the 22 GHz image with a resolution of 
1.4 X 0.95 mas: its maximum deconvolved size is 0.7 mas. GMR-A is slightly extended in 
the 15 GHz image, which has a resolution of 2.0 ± 0.60 mas. At this frequency the source is 
equally well-modeled as a circular Gaussian of deconvolved size 0.9 mas, or as two point-like 
components separated by 0.8 mas and of flux density 8.2 ± 0.3 and 3.3 ± 0.3 mJy. 

No source was detected at either 15 or 22 GHz in the first VLBA epoch at a level of 
5(7. The rms noise in the 15 and 22 GHz images is 1.3 mJy and 1.6 mJy, respectively. 
The VLA fluxes on the same day were 18.5 ± 1.7 mJy and 23.7 ± 0.6 mJy, respectively. A 
number of factors may contribute to the nondetection of GMR-A in the flrst epoch, and 
to the discrepancy in flux density in the second epoch. Instrumental errors due to phase 
decorrelation and amplitude calibration may contribute substantially. Additionally, we know 
that the source can vary significantly on a time scale of hours. Another plausible explanation 
is that some of the fiux is in a low-surface brightness "halo" which is resolved out by the 
VLBA. We discuss this possibility, which is supported by earlier VLA observations, in more 
detail in 53.1. 



2.4. Near Infrcired Photometry 

Observations of GMR-A were made on 2003 January 22 at 6 UT with the Keck I 10-m 
telescope on Mauna Kca, Hawaii, and on 2003 January 24 at 1 UT with the Cerro Tololo 
Interamerican Observatory Blanco 4-m telescope at Tololo, Chile. 

The Keck observations used the facility near-infrared camera (NIRC; Matthews & Soifer 
1994), which is equipped with a 256 x 256 pixel Santa Barbara Research Corp. InSb array, 
with a pixel size of 0.151". The images were taken at an airmass of 1.26, and under seeing 
conditions of 0.9". We fiat-fielded the data according to standard procedures, and replaced 
bad pixels with the median of surrounding pixels. The absolute calibration of the images 
was set by observing the IR standard star SJ9118 (Persson et al. 1998) at a similar airmass 
(1.28). Photometry in J, H, and Kg bands was obtained for GMR-A and two other stars 
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in the field. We simply integrated the intensity over the star and subtracted the average 
background emission. Corrections from K5 to K band were less than 0.1 mag. 

We obtained additional IR photometry from the Blanco 4-m Telescope at the Cerro 
Tololo Interamerican Observatory with the facihty wide-field camera Infrared Side Port Im- 
ager (ISPI; Probst et al. 2003). The instrument uses a Rockwell 2048x2048 pixel HgCdTe 

detector with a 10.2' field of view. Conditions were photometric with seeing of 1.6" FWHM. 
GMR-A was too bright at fC^-band to be observed in the shortest possible integration time, 
so we de-focused the telescope in order to obtain unsaturated images for this filter. We 
observed the standard star SJ 9116 from Persson et al. (1998) for photometric calibration, 
also with the telescope defocused for A'^-band. The data were reduced in a standard fash- 
ion. We subtracted an average bias frame from the images. We constructed flat fields from 
images of the lamp-illuminated interior of the dome. The individual images were registered 
via cross-correlation and stacked to form a final mosaic. Aperture photometry was used to 
measure the object and standard star fiuxes. 

Comparison between the Keck and CTIO results for GMR-A and two other field stars 
(Table 4) suggests that the photometric accuracy is on the order of 0.1 mag. 

2.5. Near Infrared Spectroscopy 

We observed GMR-A with the Keck II 10-m telescope on 2003 February 6, using the 
facility near-infrared echelle spectrometer (NIRSPEC; McLean et al. 1998). We obtained 
high-resolution {R ~ 21, 000), cross-dispersed spectra in the wavelength range 2.10-2.40 iim 
using the NIRSPEC-7 order-sorting filter. The star was imaged in two nods along the 0'.'432 
X 24" slit, separated by 12" (spatial resolution was seeing limited at ~ 1"). Each integration 
was 300 seconds, providing a total time of 600 seconds on the star and a signal-to-noise 
ratio of approximately 160. The spectra were dark subtracted, flat-fielded, and corrected for 
cosmic rays and bad pixels. The curved echelle orders were then rectified onto an orthogonal 
slit-position/wavelength grid based on a wavelength solution from sky (OH) emission lines 
and He/Ar/Ne/Kr arc lamps. Each pixel in the grid has a width of 6X = 0.031 nm. We 
subtracted sky emission by taking the difference of the two nods and then fitting third-order 
polynomials to the 2D spectra column-by-column to remove any temporal variation. 

The star's spectrum was extracted using a Gaussian weighting function matched to the 
wavelength-integrated profile of each echelle order. To correct for atmospheric absorption, 
we observed a BOV star, HD 36512; both GMR-A and the cahbration star were observed 
at an airmass of ~ 1.2. To account for photospheric absorption features and continuum 
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slope, the calibration star's spectrum was divided by a spline function fit. The resulting 
atmospheric absorption spectrum was then divided into the GMR-A spectra. 

It is difficult to separate stellar from nebular emission for GMR-A. Br 7 emission at 
2.166 ijlui fills the slit with a velocity width of 29.7 km s~^ and an equivalent width of ~ 6 A. 
Subtraction of the mean off-source spectrum from that of GMR-A reveals a weak emission 
fine with FWHM ~ 23 km s"^ and EW ~ 0.8 A. 



3. Discussion 

3.1. Radio, Infrared, and X— ray Properties 

The mm- wavelength flare source has counterparts at radio, infrared and X-ray wave- 
lengths. The object was first identified (as "source A") by Garay et al. (1987) from VLA 
observations obtained in 1981; its fiux densities were 10.0 ± 1.0 mJy and 11.0 ± 1.0 mJy 
at 5 and 15 GHz, respectively. Later VLA observations by Felli et al. (1993a) gave fiux 
densities of 5.5 ±0.5 mJy and 5.4 ±0.3 mJy at 1.4 and 15 GHz. Both studies found evidence 
that the source was slightly extended; FeUi et al. found a deconvolved size of 0.19 x 0.14 
arcsec. Prom observations at 13 epochs over an 8-month period in 1990, FeUi et al. (1993b) 
found that the fiux density of GMR-A was highly variable at 5 and 15 GHz, with a fiat 
spectrum. Our VLBA observations, in combination with the minimum fiux density in the 
VLA lightcurve, suggest that GMR-A consists of a constant, extended, fiat spectrum source 
with a fiux density ~ 5 mJy which is resolved out by the VLBA and a variable, compact 
source that is detected by the VLBA. 

As shown in Figure 6, the radio spectral index was highly variable in the weeks following 
the initial BIMA detection. On 22 January, 2 days after the fiare was detected at BIMA, the 
emission had a fiat spectrum between 22 and 43 GHz {a ^ —0.1 for S oc u"), but appeared 
to cut off sharply below 22 GHz. Two days later, as the first outburst cooled, the spectrum 
remained peaked at roughly 22 GHz, with a higher frequency spectral index of —0.5. One 
day later, on 25 January, Nobeyama observations (Nakanishi et al. 2003) detected the source 
at 98 and 147 GHz; the spectrum was inverted {a — 1.4), implying that the synchrotron 
self- absorption frequency of the plasma was greater than ~ 100 GHz. On 7 February the 
source appeared to have a fiat spectrum from 15 to 86 GHz. Much later, on 27 March, 
the spectral index is —0.7, although the absence of an 86 GHz observation makes this steep 
spectrum somewhat uncertain. 

Gircular polarization is clearly detected in many epochs (Table 3). When either the 
circular polarization or the total intensity is not detected at the 2 — a level, we report an 
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upper limit to the absolute value of the circular polarization. The mean fractional circular 
polarization (measured as V/I) is -3.4 ± 0.4%, -3.9 ± 0.2% and -6.6 ± 0.6% at 15, 22 and 
43 GHz, respectively. Nevertheless, there is clear time variability in the circular polarization: 
in particular, we did not detect circular polarization near the peak of the flare. The emission 
is never linearly polarized. We searched for linear polarization but did not detect it in any 
epoch with limits ~ 1%. 

The near infrared counterpart of GMR-A is bright, and not variable. GMR-A is po- 
sitionally coincident with source #573 in the infrared survey of Hillenbrand and Carpenter 
(2000). Our infrared photometry gives H and K magnitudes consistent with those measured 
by Hillenbrand and Carpenter (Table 4). Furthermore, the intensive survey for variability in 
the Orion nebula by Carpenter, Hillenbrand & Skrutskie (2001) did not detect this object 
as significantly variable. 

In contrast, the X-ray counterpart of GMR-A is characterized as a flare source by 
Feigelson et al. (2002); it is substantially variable on a timescale of less than 12 hours as 
well as over six months. The spectrum is consistent with an intrinsic X-ray luminosity 
— 10^^-^ erg s~^ attenuated by a gas column density Nh — 10^^-^ cm"^. This luminosity 
ranks it among the brightest 10% of the X-ray sources in the Orion nebula. Getman et al. 
(2003) reported a signiflcant X-ray flare (xlO) at its position. Their preliminary light curve 
shows the flare beginning ~ 2 days before the detected millimeter wave flare, and continued 
through the time of the radio flare. Of course, we arc not able to accurately determine 
the onset of the flare at millimeter wavelengths because of the scarcity of observations. The 
X-ray flux was declining at the time of the millimeter flare. Modeling of the X-ray spectrum 
by Getman et al. during the flare indicates an absorption column density comparable to that 
found earlier. 



3.2. Vciriability Timescale of the Radio Emission 

To study the flux variability during the onset of the radio flare, we broke up the 8-hour 
long visibility data obtained by BIMA on the night of 2003 January 20 into 1-hour blocks. 
The flare brightened by a factor of ~ 4 to 8 during the track (Figure 2), while the fluxes of 
IRc2 and BN remained stable. The flux density of the flare source was fit in each interval 
after subtracting the contributions from point sources at the positions of IRc2 and BN with 
flux densities of 45 and 68 mJy respectively. Because the BIMA receivers are sensitive to 
linear polarization in the vertical direction, a concern is that the observed variability may 
not be intrinsic, but introduced by linearly polarized emission combined with changes in 
the parallactic angle during the night. The parallactic angle changed from —38° to +45° 
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during the observations. If this were the cause of the variabihty, the emission would have 
to be > 45% hnearly polarized, which is highly unusual. A broad range of position angles 
and polarization fractions were used to model the light curve, obtaining a minimum reduced 
~ 17, suggesting that linearly polarized emission cannot reproduce the observations. 
Furthermore, linear polarization is undetected in all subsequent VLA observations. It is 
very unlikely that the 86 GHz flux is initially strongly polarized while two days later the 43 
GHz flux is not. This fact, together with the poor quality of the fit, strongly argue that the 
observed 86 GHz lightcurve is due to intrinsic variability and not to polarization. 

We searched for short-term variations in the VLA data, as well. Most of the observations 
were split into two or more blocks probing timescales that range from 0.5 to 1.5 hour. There 
was no evidence for significant variability on these timescales. 

3.3. Identification of GMR-A as a Young Stellar Object 

To determine the spectral type of GMR-A, wc compared its near-IR spectrum with 
the medium-resolution fT-band spectral atlas of Wallace & Hinkle (1997). The relative 
strengths of the Al and Mg lines near 2.11 /xm, as well as the relative strengths of the Ti 
and Si lines near 2.18 jsm, indicate a spectral type of K5. The shallow CO bandheads are 
consistent with a dwarf star. On this basis, we identify GMR-A as a spectral type K5V. 

The ^^CO (v=0 2) rovibrational band (extending redward from the bandhead at 
2.29 ixm) allows us to measure the radial and rotational velocities of the star, visr and 
vsini, respectively. As a template for our rotational model, we used a high-resolution FTS 
sunspot umbral spectrum (Wallace & Hinkle 2001) with Tg// — 4100K. We convolved the 
sunspot spectrum with a rotational broadening profile from Gray (1976), and convolved the 
result with the measured profile of an OH line observed with NIRSPEC to account for the 
instrumental resolution. Iterating over a range of rotational velocities, we find v smi = 23±2 
km s^^ and Visr = —4 ± 5 km s^^. The quoted error in vsini reflects the broad nature of 
the minimum in the residuals for models near 23 km s~^. Our determination of Vi^j. suggests 
that GMR-A is not embedded in the molecular cloud, for which the radial velocity from NH3 
measurements is visr = 9 ± 1 km s"-*^ (Batrla et al. 1983), but rather lies behind it. 

Despite the difficulties in subtracting the background emission from the nebula, the 
observed Br 7 emission is clearly substantially weaker and narrower than is typical in classical 
T Tauri stars (Folha & Emerson 2001), indicating that this is probably a weak hne T Tauri 
star (WTTS). This designation is also consistent with strong radio variability, as WTTS are 
the most luminous and most variable of radio YSOs. 
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Taken together, the X-ray data, IR photometry, and IR spectroscopy point to a clear 
identification of GMR-A as a YSO inside, or perhaps behind, a molecular cloud. However, 
the identification of GMR-A as a K5V is not consistent with the high luminosity implied 
by the very bright K magnitude and the large extinction. We discuss here the different 
scenarios that might resolve this, including different degrees of extinction, the presence of 
an infrared excess, and classification as a K5III star, but none of these are fully satisfactory. 
In the simplest approach, we use the X-ray absorption column depth and Ry — 3.1 to find 
the extinction Ay — 21.4 which implies — 2.4 mag (Rieke & Lebofsky 1985; Binney & 
Merrifield 1998). Using color and bolometric corrections appropriate for K5V stars, we find 
L ~ 30^0. But this result is almost certainly wrong for two reasons: (1), this luminosity 
exceeds the luminosity of the stellar birth line, which is < IOL0 at this effective temperature 
(Palla & Stabler 1999); and (2), the photospheric H-K color of 1.0 after applying reddening 
correction is not appropriate for a K5V star. Assuming a photospheric color H-K=0.1, which 
is the case for K5V stars, we find Ay = 35 and Ak = 3.9. A similar extinction is obtained 
with examination of the J-H and H-K colors (e.g., Haisch, Lada & Lada 2000). While this 
satisfies the spectroscopic constraints, it produces a luminosity of 100 Lq, again in excess 
of the stellar birth line. This inconsistency suggests an excess of infrared emission possibly 
due to a disk. Assuming the reddening calculated by the X-ray absorption column depth 
(Ak = 2.4) and a photospheric color H-K=0.1, then we find an H-K excess of 0.9 mag and 
a K-band excess of 1.7 mag (Hillenbrand & Carpenter 2000). This reduces the bolometric 
luminosity to 6 Lq. While this satisfies the stellar birth line luminosity limit and the colors 
for these stellar classes, the presence of a strong infrared excess is more consistent with a 
classical T Tauri star than with a WTTS, although some WTTS do have an IR excess (Strom 
et al. 1989). Considering the same scenarios for K5 giant or supergiant stars leads to the 
same problems. 

Could GMR-A be something other than a T Tauri star? A possibility is that the flaring 
source active in radio and X-rays is not the object dominating the IR luminosity, hence the 
lack of infrared activity. There is some precedent for faint companion stars producing intense 
variable radio emission: astrometric radio and optical measurements of 9^ Ori A have shown 
that the object responsible for radio flaring is actually part of a triple system (Garrington 
et al. 2002). The "active" companion, however, would have to be several magnitudes fainter 
than the main K5V star even during its flaring state, to account for its non detection in the 
IR spectrum. Adaptive optics observations in the infrared would be necessary to identify any 
faint companions (e.g., Duchene et al. 2003). Other potential radio transient source types 
are brown dwarfs and late-type main sequence stars (e.g., Berger 2002). These sources, 
however, are too faint by many orders of magnitude to account for the quiescent or flaring 
radio emission from GMR-A. Evolved stars such as RS Cvn binaries are unhkely to be 
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associated with the Orion complex. In the following section we show that the magnetic field 
measurements also support the conclusion that GMR-A is a T Tauri object. 

3.4. Zeeman Measurement of Magnetic Field Strength 

We used the echelle infrared spectrum to search for Zeeman broadening of lines in 
the stellar photosphere, which would indicate the presence of a strong magnetic field. We 
selected the strong Til fines at A = 2.17886 iim and A = 2.19033 iim present in echelle 

order 35 in our spectra to carry out this study. These transitions have spectroscopic terms 
(3d^ 5p^ ^4p^ i) and (3d^ ^4p^ ^Ds) respectively, resulting in very similar effective 
Lande (yf-factors, 1.667 and 1.833. Because the factors of the upper and lower energy lev- 
els in either transition are not precisely equal, these lines formally exhibit the "anomalous" 
Zeeman effect and they are broken up into seven components by the magnetic field. In prac- 
tice, since the Lande factors are almost equal, the emission is separated into approximately 
3 equal components. 

We modeled the observed emission using 3 spectral lines with the shape of the intrinsic 
instrumental profile convolved with the rotational profile corresponding to the measured 
value of vsini — 23 km s~^. Thus, we assume that the only sources of line broadening 
are instrumental, rotational, and Zeeman, an assumption that seems justified in view of the 
large vsini. These models were computed for a grid of Zeeman splittings, AA, and fine 
amplitudes. We determined the best model by evaluating the reduced oi the residuals. 
The error was estimated by finding the values of the parameters that increased the reduced 
by one. Figure 7 shows the results of the fit for the A = 2.17886/im transition. The measured 
Zeeman splittings are AA = 0.14 ±0.07 nm (x^ = 1.4) and AA = 0.09 ±0.05 nm (x^ = 1.7), 
corresponding to magnetic fields 3.8 ±2 kG and 2.2 ±1.1 kG. Thus, our best estimate of the 
photospheric magnetic field traced by the Til lines is 2.6 ± 1.0 kG. Observations of T Tauri 
objects find similar magnetic field strengths (e.g., Johns-KruU, Valenti & Koresko 1999). 

3.5. Radio Emission Mechanism 

Several lines of evidence indicate that the radio emission mechanism is nonthermal: the 
compactness of the emission detected using the VLBA, the presence of circular polarization, 
the strongly variable spectral index, and the rapid rise and decay time of the total intensity. 
We can set a lower limit to the brightness temperature based on the fiux measured by the 
VLBA observations. We find T,, > 5 x 10^ K and 7], > 3 x 10^ K at 15 and 22 GHz. Electrons 
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in equilibrium with this temperature will have velocity > 0.07c, or an electron Lorentz factor 
7e > 1.003. 

The presence of high fractional circular polarization rules out thermal bremsstrahlung 
and synchrotron radiation as emission mechanisms, and indicates that the dominant emission 
is cyclotron radiation generated by mildly relativistic electrons in the presence of strong 
magnetic fields (e.g., Giidel 2002). The observed mean fractional circular polarization in the 
VLA bands is ~ —3 to —7%, with peak values greater than —10%. The absence of linear 
polarization also argues against synchrotron emission, although very tangled magnetic fields 
or a high rotation measure due to an intervening dense, magnetized plasma could depolarize 
a synchrotron source. 

The radio emission probably originates in coronal magnetic fields, which are an order of 
magnitude weaker than the photospheric magnetic field strength determined in the previous 
section. Assuming a field strength for the radio emission of 100 G, the corresponding electron 
gyro-frequency is 0.57~^ GHz. The lack of a very steep fall-off at high frequencies, char- 
acteristic of a single temperature thermal cyclotron emitter, leads us to conclude that the 
emission originates from a power-law distribution of electrons with low energy. A power-law 
distribution of energies with index S — 1.5 yields an optically-thin spectrum with a — —0.1 
above a peak frequency of 22 GHz. The characteristic size of the emitting region for the 
assumed magnetic field is 2 x 10^^ cm. For comparison, models indicate that the radius of 
a YSO with L = 3.5Lq and T^ff = 4100 K is ~ 4 x 10^^ cm, suggesting that the cyclotron 
emission is a global phenomenon for the star. 

A significant anomaly in this picture is the absence of circular polarization in the first 
VLA observation on January 22, 2 days after the millimeter wavelength flare. A possible 
mechanism to initially supress circular polarization is to invoke an initial outburst of relativis- 
tic electrons, leading to synchrotron rather than cyclotron emission. Rapid cooling through 
expansion or X-ray radiation would then drive the radio source into the cyclotron regime. In 
this scheme one expects that the early radiation would be linearly polarized. We have seen 
that VLA observations on day 22 failed to show any hnear polarization. As discussed above, 
however, there are mechanisms that could significantly depolarize a synchrotron source. An 
alternative explanation for the initial lack of circular polarization assumes the existence of 
two oppositely polarized sources associated with the star. Imaging of some flaring radio 
stars has shown such oppositely polarized lobes (Mutel et al. 1998). Circular polarization 
disappears when these lobes are equally balanced and appears when one lobe comes to dom- 
inate the emission. Without more detailed modeling and/or imaging of the source we cannot 
discriminate between these two possibilities. 

The source reflared several times during the course of the follow-up observations, albeit 
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never to the level of the discovery data. In fact, the light curve appears to consist of multiple 
flares spread over ~ 70 days, with individual flares having a rise and decay timescale of 
days or shorter. During the BIMA discovery, we measure a rise time of ~ 1 hour. Since this 
is the duration of most of the VLA observations, we cannot clearly discern whether there is a 
similar variability timescale in those data. The timescales for radiation losses for a cyclotron 
source are proportional to B~^, and are ~ 1 day for B = 100 G, which is comparable to 
what is observed. 



The luminosity of this flare is comparable to that of the brightest stellar radio flare 
ever detected. The peak luminosity at 86 GHz during the flare was 4 x 10^^ erg s^^Hz^^, 
only 50% less than the peak flare of the current record-holder, the FK Gom-type giant 
star HD32918 (Slee et al. 1987, Bunton et al. 1989). The luminosity of the GMR-A flare, 
however, was more than an order of magnitude higher than detected for any YSO. Typical 
maximum radio luminosities for YSOs are ~ 10^^ erg s^^Hz"^ (Giidel 2002). Since GMR-A 
flared at millimeter wavelengths, its total luminosity was L,^ ~ 3 x 10^*^ erg s"-*^, about two 
orders of magnitude higher than that of centimeter wavelength sources. 

The ratios of X-ray and radio emission for GMR-A during both its quiescent and flaring 
states approximately follow the correlation found for active stars and solar flares (Giidel 2002) 



This correlation was derived from radio observations at centimeter wavelengths (typically 
8 GHz): the precise value of the coefficient is unknown for measurements at millimeter 
wavelengths. To compute the luminosities, we adopted the quiescent, absorption-corrected 
X-ray luminosity — 10^^'^ erg s~^ measured by Chandra. During the flaring state the 
Chandra observations show the count rate increasing by a factor ~ 10, leading to a prediction 
of 5 X 10^^^^ erg s^^Hz"^, very similar to that observed at 86 GHz. These data support 
the conclusion that the flare is caused by magnetic activity common to other radio stars. 



We can parametrize the rate of millimeter wavelength flares in the Orion nebula as 



3.6. The Radio to X— ray Luminosity Correlation 



(1) 



3.7. Rate of Flaring Activity in the Orion Nebula 
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Nonthermal activity often follows a power-law index a — 1. This is the first fiare brighter 
than 100 mJy detected by BIMA in ~ 40 observations of the Orion nebula, setting A = 
0.01 — 0.1 day^^. This value of A is consistent with the 22 and 43 GHz observations described 
in this paper, as well as the 15 GHz dataset of Fclli et al. (1993b). Both of these datasets 
show that GMR-A has a flux density of 10 mJy about 10 times as frequently as it has a 
flux density ~ 100 mJy. Following Equation 2, observations of Orion at a 0.1 mJy detection 
threshhold will find ~ 10 — 100 fiares from YSOs per observation. In fact, the monitoring 
campaign of FeUi et al. (1993b), with a sensitivity ranging between 0.2 and 1.2 mJy, identified 
13 variable sources. The proposed Atacama Large MiUimeter Array will achieve a sensitivity 
of 0.1 mJy or better in a 1 minute integration (http://www.alma.nrao.edu). Thus, a short 
ALMA observation with a sensitivity of 10 /uJy may find ~ 100 — 1000 fiares from YSOs, 
which is an appreciable fraction of the total number of objects in the region. 

The same argument can be made using the X-ray luminosity function. With a radio 
sensitivity of 3 x 10^^ erg s~^ Hz~^ (equivalent to lO/xJy), the L^/Lr correlation predicts that 
we will detect all stars with L^. > 10^°'^^"^. This is consistent with the mean logL^- = 29.4±0.7 
of the Chandra sample, suggesting that we would detect half of all X-ray objects in the Orion 
nebula (Feigelson et al. 2002). Approximately half of the 1075 objects in the Chandra sample 
are identified as variable on short and/or long timescales, indicating again that ALMA can 
expect to find hundreds of variable objects in a few hours. 

Variable sources violate the assumption of a constant sky that underlies synthesis imag- 
ing. The potentially large number of variable sources in the Orion nebula and other star- 
forming regions introduces a complication not yet considered for imaging with ALMA. If 
not adequately accounted for, these variable sources will introduce a dynamic range limit to 
images. Consider a source with a fiare of fiux density S^ax over the quiescent fiux density 
of Sq that lasts for tmax « tabs, whcrc tabs is the total length of the observation. An ar- 
ray will have a peak fractional sidelobe level rjmax that is a function of array configuration, 
source position and tmax- To first order, this will introduce an imaging error f]maxSmax after 
deconvolution (e.g., using CLEAN). For ALMA, rj ~ 2% for observations of length > 1 hour 
(e.g., Boone 2002). Thus, a brief 1 mJy flare will introduce an error > 20 fiJj in an image 
that could have a theoretical rms of 10 /xJy Multiple sources and multiple flares complicate 
this picture but are unlikely to reduce the error. The solution, in principle, is to image and 
deconvolve on a timescale comparable to the flaring timescale. 
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4. Conclusions 

We have described here the serendipitous discovery and subsequent follow-up observa- 
tions of a flaring millimeter wavelength source in the Orion star-forming cluster. The source 
appears to be a weak-line T Tauri star highly obscured by dusty material. Although pre- 
viously known as a variable cm radio source (GMR-A), the magnitude of the flare and the 
measurement of its spectrum to millimeter wavelengths are novel. This is one of the most 
luminous radio flares observed, and the most luminous flare from a YSO. Future observations 
with greater sensitivity can be expected to detect many more of these objects. 

The discovery of a compact source, detectable by the VLBA, in the Orion nebula raises 
the possibility of the measurement of trigonometric parallax to this star-forming region. It 
would be the most accurate distance measurement for Orion, achieving a parallactic error 
~ 100 //as equivalent to a distance error of ~ 2% with only a few measurements. If the 
quiescent source is detectable, or if the source flares again, then two more observations 
are sufficient to separate the proper motion and parallax of this source. Monitoring and 
detection of other transient sources in this region and other galactic star-forming regions 
could significantly reduce the uncertainty in key star-forming parameters. A statistical 
approach is possible through VLBI astrometry of the > 10 other nonthermal sources in 
Orion previously identified by FeUi et al. (1993a). Finally, this may allow a 3-dimensional 
probe of the positions and velocities in the star forming cluster, giving vital clues to its 
history. 
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Fig. 1. — BIMA 86 GHz continuum image of Orion-KL taken on 2003 January 20, showing 
the flare source at upper right. The synthesized beam is 0'.'93 x 0'.'56 at P.A. 43°. The halftone 
scale ranges from 12 to 60 mJy beam~^. The map has not been corrected for attenuation 
by the 2^2 FWHM primary beam, shown by the circle. Source I, the radio source associated 
with IRc2, is at the center the map. The faint arc just below source I is dust emission from 
the Orion hot core. The Becklin-Neugebauer object is 10" NW of source I. 
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Fig. 2. — BIMA flux densities at 86 GHz with a model of linearly polarized flux assuming 
p — 50% and a position angle of 107 degrees. While the polarization model follows the trend 
of the data, it is not an adequate fit for any position angle and polarization fraction. 
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Fig. 3. — Light curve at millimeter, radio and X-ray wavelengths showing early evolution 
of the flare along with dates of infrared photometry ("IR" arrows), VLBA observations 
("VLBA" arrows) and infrared spectroscopy ("IR Spec" arrow). Symbols are as given in 
the legend. The downward-facing red triangles ( "BIMA 86 GHz UL" ) indicate 3 — a upper 
limits. The X-ray light curve is from the description in Getman et al. (2003) and is given 
in arbitrary flux units. 
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Fig. 4. — Full light curve at millimeter and radio wavelengths showing the long term evo- 
lution. Symbols are as given in the legend. The downward-facing red triangles ("BIMA 86 
GHz UL" ) indicate 3 — 0" upper limits. 
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Fig. 5 — VLBA images of GMR-A at 22 GHz obtained on 2003 January 29. Contours are -3, 
3, 6, 12, and 24 times the rms noise levels of 0.70 mJy beam~^ at 22 GHz. The synthesized 
beam is shown in the lower left hand corner of the image. 



-24- 
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Fig. 6. — Spectra of GMR-A obtained on 2003 January 22 (squares, heavy solid line), January 
24 (circles, dotted line), January 25 (crosses, dashed line), February 17 (stars, dot-dashed 
hne) and March 29 (diamonds, hght sohd hne). Upper hmits are 3 — cr. 1 — a errors are 
plotted where they exceed the symbol size. The 97 and 148 GHz data of 25 January 2003 
are from Nakanishi et al. (2003). 
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Fig. 7. — Three component Zeeman fit of tfie 2.17886 //m transition of Ti I in GMR-A. The 
three Zeeman components, broadened by the instrumental profile and convolved with the 
rotational profile for v sini = 23 km s^^, are shown by dashed lines. A nearby Si I transition 
has been modeled in the same manner and included in the fit (dotted line). The sum of all 
the components is shown by the thick gray line, superimposed on the data. A high resolution 
spectrum of a sunspot with effective temperature similar to a K4 star and a strongly split 
Ti I fine is shown for comparison. The splitting measured by the fit is AA = 0.14 ± 0.07 nm, 
corresponding to a magnetic field 3.8 ± 2.0 kG. 
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Table 1. Position of GMR-A 



Observation 


Observing Band 


RA 


Dec. 






(J2000) 


(J2000) 


VLBA 


15 GHz 


05 35 11.802691 ± 0.000001 


-05 21 49.24660 ± 0.00004 


VLBA 


22 GHz 


05 35 11.802695 ± 0.000001 


-05 21 49.24660 ± 0.00003 


VLA 


15, 22, 43 GHz 


05 35 11.801 ± 0.001 


-05 21 49.27 ± 0.02 


BIMA 


86 GHz 


05 35 11.80 ± 0.006 


-05 21 49.2 ± 0.1 


IR 


2.2 //m 


05 35 11.81 


-05 21 49.3 


X-Ray 


1— 10 keV 


05 35 11.8 


-05 21 49 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation 
are degrees, arcminutes, and arcseconds. 
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Table 2. BIMA Fluxes for GMR-A 



Day 5*86 
(mJy) 



-8.29 


< 6 


11.29 


< 12 


14.29 


< 15 


20.12 


42 ± 9 


20.17 


39 ± 9 


20.21 


18 ± 7 


20.25 


122 ± 7 


20.29 


136 ± 7 


20.33 


146 ± 7 


20.38 


162 ± 8 


23.23 


< 25 


24.23 


< 13 


30.38 


< 27 


35.23 


< 6 


38.16 


11 ± 3 



Note. — 



Upper limits are 3a corrected for attenuation by the primary beam. 
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Table 3. VLA Fluxes for GMR-A 



UT Day ^8.4 


'S'ls 


'S'22 


'S'43 


^15 


V22 


F43 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 



22.23 


< 33 








85.7 




0.7 


76.9 


± 


1.0 






< 


1.4 


< 


2.8 


24.23 




18.5 




1.7 


23.7 




0.6 


16.6 




0.5 


< 


2.6 


-1.2 ± 


0.3 


-1.9 ± 


0.5 


27.26 




43.6 


± 


1.6 


43.7 


± 


0.5 


35.7 


± 


0.4 


< 


1.4 


-2.2 ± 


0.2 


-3.4 ± 


0.5 


30.24 




32.4 


± 


1.7 


30.1 


± 


0.8 


26.6 


± 


0.5 


< 


2.0 


< 


1.0 


< 


2.3 


31.25 




33.8 


± 


1.6 


34.7 


± 


0.6 


25.9 


± 


0.5 


< 


1.9 


-0.6 ± 


0.2 


< 


1.9 


32.17 




27.3 


± 


1.6 


19.2 


± 


0.6 


15.3 


± 


0.4 


< 


1.8 


-1.1 ± 


0.2 


< 


2.1 


33.21 




41.6 


± 


1.6 


46.8 


± 


0.6 


28.0 


± 


0.5 


-1.3 ± 


0.3 


-1.8 ± 


0.2 


-1.9 ± 


0.4 


34.21 


< 31 


17.8 


± 


1.0 


19.9 


± 


0.3 


10.9 


± 


0.3 


-1.2 ± 


0.2 


-1.0 ± 


0.1 


-0.9 ± 


0.3 


36.23 




10.1 


± 


2.5 


8.7 


± 


1.3 


8.6 


± 


0.5 


< 


1.8 


< 


3.1 


< 


1.9 


37.98 




9.9 


± 


3.2 


10.2 


± 


1.1 


5.5 


± 


1.1 


< 


4.0 


< 


1.2 


< 


1.6 


42.13 




17.7 


± 


2.1 


7.1 


± 


0.6 




< 


1.0 


< 


1.4 


-0.6 ± 


0.1 


< 


1.3 


45.04 


< 36 


7.5 


± 


2.1 




< 


2.9 


2.6 


± 


0.5 


< 


1.8 


< 


1.8 


< 


1.6 


48.02 


< 33 


27.4 


± 


1.9 


6.6 


± 


0.7 




< 


1.7 


< 


3.5 


< 


1.5 


< 


1.3 


53.17 


< 21 


10.2 


± 


1.6 


9.4 


± 


0.6 


9.2 


± 


0.3 


-1.2 ± 


0.2 


< 


0.9 


-1.0 ± 


0.2 


55.10 




7.3 


± 


1.8 


4.7 


± 


0.7 


6.7 


± 


0.3 


< 


1.5 


< 


0.9 


< 


1.2 


59.12 


< 26 


25.8 


± 


1.9 


10.2 


± 


0.8 




< 


1.5 


-0.7 ± 


0.2 


-1.0 ± 


0.1 


< 


1.1 


76.06 


< 18 


24.9 


± 


1.1 


19.5 


± 


0.5 


15.9 


± 


0.3 


< 


1.2 


-0.5 ± 


0.1 


-1.9 ± 


0.3 


83.04 


< 21 


10.9 


± 


1.0 


9.1 


± 


0.4 


5.0 


± 


0.2 


-0.9 ± 


0.1 


-0.9 ± 


0.1 


< 


1.2 


88.96 


< 36 


8.8 


± 


1.4 


6.1 


± 


0.5 


3.2 


± 


0.4 


< 


1.0 


-0.9 ± 


0.1 


-0.6 ± 


0.1 



Note. — Sf and Vf are the total intensity and circularly polarized flux density at frequency 
/. Upper limits given for V} are 2(7 limits on the absolute value. UT days are given for the 
year 2003. 
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Table 4. Infrared Photometry 



Star 


Date 


Telescope 


J 


H 


K5 


Reference 


GMR-A 


02/09/99 


Keck/NIRC 




11.9 


9.6 


HC2000 #573 


GMR-A 


01/22/03 


Keck/NIRC 


16.0 


11.98 


9.61 


This paper 


GMR-A 


01/24/03 


CTIO 


15.8 


11.9 


9.7 


This paper 


Star W 


02/09/99 


Keck/NIRC 




10.5 


9.9 


HC2000 #554 


Star W 


01/22/03 


Keck/NIRC 


11.1 


10.38 


9.75 


This paper 


Star W 


01/24/03 


CTIO 


11.2 


10.3 


9.8 


This paper 


Star S 


02/09/99 


Keck/NIRC 




12.0 


11.6 


HC2000 #555 


Star S 


01/22/03 


Keck/NIRC 


12.7 


12.1 


11.5 


This paper 


Star S 


01/24/03 


CTIO 


12.7 


12.0 


11.6 


This paper 



Note. — Star S is the bright K band object 6" South of GMR-A. Star W is a bright K 
band object 6" South and 20" West of GMR-A. HC2000 refers to Hillenbrand & Carpenter 
(2000). 



